Thermal rounding exponent of the depinning transition of an elastic string in a 

random medium 
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We study numerically thermal effects at the depinning transition of an elastic string driven in a 
two-dimensional uncorrelated disorder potential. The velocity of the string exactly at the sample 
critical force is shown to behave as V ~ with tp the thermal rounding exponent. We show 
that the computed value of the thermal rounding exponent, ip = 0.15, is robust and accounts for 
the different scaling properties of several observables both in the steady-state and in the transient 
relaxation to the steady-state. In particular, we show the compatibility of the thermal rounding 
exponent with the scaling properties of the steady-state structure factor, the universal short-time 
dynamics of the transient velocity at the sample critical force, and the velocity scaling function 
describing the joint dependence of the steady-state velocity on the external drive and temperature. 

PACS numbers: 64.60.Ht, 75.60.Ch, 05.70.Ln 



I. INTRODUCTION 

The understanding of the static and dynamic proper- 
ties of elastic interfaces in disordered media has direct 
impact on different fields on condensed matter physics. 
Among a large variety of systems one can mention mag- 
netic ]l|-|I| or ferroelectric [jj Q domain walls, contact 
lines 0|) fractures @, vortex lattices fiol - fl^ . charge 
density waves [l3j], and Wigner crystals [14| . as paradig- 
matic examples. Since the effect of the disordered media 
in all these systems is non-trivial, an important question 
is how these elastic objects respond to an external drive. 

When the temperature is zero, there exists a critical 
force value F c such that the steady state velocity of the 
center of mass of the interface is zero below F c and is 
finite above it. This is due to the complex interplay be- 
tween disorder and external force: the interface accom- 
modates within the disorder energy landscape and a finite 
energy barrier must be overcome by the external force in 
order to generate a net movement. Therefore a finite 
force value has to be set to have an infinitcsimally small 
finite velocity. This is the so called depinning transition. 
If the critical force value is approached from above, the 
velocity vanishes as V ~ (F — F C )P for a thermodynamic 
system, with j3 the depinning exponent. Concomitant 
with the power-law decrease of the velocity is the di- 
vergence of a characteristic length as £ ~ (F — F c )~ u , 
with v the correlation length exponent. This depinning 
correlation length gives the typical size of the correlated 
displacement (or avalanche) that makes the interface ad- 
vance in the direction of the external force. The finite 
force threshold, the critical decrease of the velocity order 
parameter and the divergence of the typical length scale 
led to propose a description of the depinning transition 
using tools from standard critical phenomena [l5[ . More 
recently however, the analysis of the low-temperature 
averaged steady-state geometry has shown that no di- 
vergent steady-state correlation length-scale exists ap- 
proaching the critical force from below, thus breaking the 



naive analogy with standard phase transitions, where two 
divergent length-scales are expected above and below the 
critical point [HI, [l?} • 

When the temperature is finite there is no sharp tran- 
sition between zero and finite velocity regimes. Even at 
forces much smaller that the critical value the interface is 
able to move since thermal activation is enough to over- 
come the effective energy barriers generated by the disor- 
der. This regime, F <C F c , is the creep regime, and it is 
characterized by a stretched exponential dependence of 
the velocity with the inverse of the external force (l8l - [23| . 
On the other hand, at forces around the critical value, 
F « F c , a finite temperature value smears out the tran- 
sition, which is no longer abrupt. This thermal rounding 
of the depinning transition can be characterized, exactly 
at the critical force F — F c , by a power- law vanishing of 
the velocity with the temperature as V ~ , with ip the 
thermal rounding exponent (23430| . 

The values of the different exponents characterizing 
the depinning transition are universal in the sense that 
their values depend on few parameters of the system such 
as the range of the intrinsic elasticity, the dimensional- 
ity of the problem, and the correlated structure of the 
disorder. For the experimentally relevant case of 1 + 1 
dimensional elastic interfaces moving in a random-bond 
disorder environment with short-range correlations and 
short-range elasticity, we have recently reported the value 
ip = 0.15±0.01 using Langevin dynamics numerical sim- 
ulations [13] ■ This value compares well with the value 
ip = 0.16 reported in Ref. [25| based in numerical simu- 
lations. However, these values are smaller than the value 
ip = 0.24 obtained using an artificial extremal activated 
dynamics [28j , which might indeed be in a different uni- 
versality class. The value ip = 0.2 was obtained using 
numerical simulations of domain wall motion with the 
random-field Ising model [2(| [2?J ■ Although it is ex- 
pected that for T > and around the depinning tran- 
sition the characteristic exponents do not depend on the 
random-bond or random-field character of the disorder, 
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this slightly larger value might be possibly ascribed to the 
anharmonic corrections to the elasticity present in the 
random-field Ising model. On the other hand, functional 
renormalization group equations at the depinning (23j al- 
low in principle to extract the thermal rounding expo- 
nents. However, in practice there are, up to now, no an- 
alytical estimates of tp, unlike the other critical exponents 
which have been computed using functional renormaliza- 
tion group up to two loops [3 if . The very existence of 
a thermal rounding, obeying a power law scaling is not 
rigorously proven, and there are indeed some models of 
depinning which exhibit at finite temperature a totally 
different type of thermal rounding [32| . It is thus crucial, 
given the uncertainty on the very type of thermal round- 
ing and certainly on the value of the thermal rounding 
exponent, to develop new methods to determine ip, and 
to check the robustness, consistency, and expected uni- 
versality of the phenomenological scaling arguments. 

Experimentally, access to the full force range relevant 
to the depinning transition has been reported in ultra- 
thin ferromagnetic layers [1, [34|. In this case, the 
thermal rounding of the depinning transition is gener- 
ated through an effective temperature dependence con- 
trolled by the relative disorder intensity among different 
samples. Indeed, it has been shown that thermal effects 
on the velocity-force characteristics can be well described 
using the value ip = 0.15 ± 0.10 

The aim of the present work is to give further numer- 
ical support to the reported value ip — 0.15, by checking 
the robustness and consistency of the scaling arguments 
applied to different observables. To this end, we show 
how this value allows to describe different measures char- 
acterizing the critical behavior of the depinning transi- 
tion: an analysis of the finite temperature structure fac- 
tor, a short-time dynamics analysis, and the analysis of 
the scaling function describing the velocity dependence 
on force and temperature around depinning for different 
disorder intensities. 



II. MODEL SYSTEM AND NUMERICAL 
SIMULATIONS 

In order to model the dynamics of one-dimensional in- 
terfaces in disordered media we use a short-range elastic 
string, as described in the following. The string is defined 
by a single valued function u(z,t), giving its transverse 
position u in the z axis. The time evolution of the string 
is given by the overdamped equation of motion 

jd t u(z,t) = cd%u(z,t)+F p (u,z)+F + T)(z,t), (1) 

where 7 is the friction coefficient and c the elastic con- 
stant. The pinning force comes from the derivative of the 
random-bond pinning potential U(u, z), i.e. F p (u, z) = 
—d u ll(u,z), whose sample to sample fluctuations are 
given by 

[U(u, z) - U(u', z')} 2 = 5{z - z') R 2 {u - u 1 ) (2) 



where R(u) stands for a correlator of range r/ [22| . and 
the overline indicates average over disorder realizations. 
Thermal fluctuations are included through the thermal 
noise term f]{z,t) which satisfies 

(r ) (z,t)) = 0, 
(r)(z, t)7](z', t')) = 2 7 T<S(i - t')6(u - «'), (3) 

where T is the temperature (with Boltzmann constant 
set to unity, fc^ = 1) and the angular brackets denote 
a thermal average. Finally, the force F on Eq. ([1} cor- 
responds to a uniform and constant external field which 
drives the string in the u direction. 

The evolution Eq. (JXJ) is numerically solved. The z di- 
rection is discretized in L segments of size Sz = 1, i.e. 
z — s- j = 0, ...,L — 1, while keeping Uj(t) as a continuous 
variable. This sets the longitudinal finite system size L. 
The equation is integrated using the Euler method with a 
time step St = 0.01. The pinning potential is modeled by 
performing a cubic spline passing through M regul arly 
spaced uncorrelated Gaussian numbers points [35l . l36| . 
which sets the transverse finite system size M. Numer- 
ical simulations are performed using periodic boundary 
conditions in both directions and using the parameters 
7 = 1, c = 1, and rt = 1. The strength of the disorder 
is given by i?o = R(0). For each disorder realization, 
i.e. for each finite size sample, the critical force F c can 
be accurately obtained using an exact algorithm, which 
also gives the critical pinned configuration of the string 
u c (z) [35[. The results presented in the following sections 
were obtained by typically averaging over 100 disorder 
configurations; the error bars being typically of the order 
of the size of the data points. 

III. VELOCITY-FORCE CHARACTERISTICS: 
SCALED VARIABLES 

In this section we will present the general features 
of the velocity-force characteristics, allowing us to de- 
fine the critical region and the scaled variables that will 
be used throughout the rest of this work. Figure Ufa) 
shows typical velocity-force curves at finite temperature, 
as obtained with the present model for L = 1024 and 
M = 5792 « Lt*», with C dep = 1.25 [13, HI] the depin- 
ning roughness exponent (see below) . Given a fixed force 
F, the velocity is computed in the steady state, which 
is typically reached within one sweep over the transverse 
size M (as detailed below, the transverse size M will be 
varied following a scaling relation with the string length 
L). Then, of the order of five sweeps over M are used to 
compute the velocity 

V=(d t u(z,t)). (4) 

The thermal average is taken by computing 200 values of 
the velocity with independent thermal noise realizations 
within this steady state regime. Different curves corre- 
spond to the same single disorder realization with inten- 
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(a) 




FIG. 1. (Color online) (a) Velocity-force characteristics for a 
single disorder realization of intensity Ro — 2 and different 
temperatures, (b) Velocity-force characteristic for tempera- 
ture T = 0.3 and different disorder intensities as indicated. 
The inset shows the scaled data according to v = V/ (mF c ) 
and / = (F — F c )/F c , where F c is the sample dependent 
critical force. Furthermore, this rescaling strongly reduces 
sample-to-sample fluctuations. 



sity Ro — 2 and increasing temperature. The characteris- 
tic critical force is indicated in the key. The lower curve, 
corresponding to T = 0, clearly presents the typical 
abrupt depinning transition: the velocity is strictly zero 
for F < F c , while it increases as (F — F c )@ for F > F c , 
where /3 < 1 is the velocity exponent. As observed, by 
increasing the temperature the T = sharp transition 
is smeared out. Although at very small temperatures 
the curves still present the curvature corresponding to 
(F — F C )P , at higher temperatures there are no clear sig- 
natures of the underlying T — depinning transition. 
Finally, at very high temperatures, when the thermal en- 
ergy is larger than the typical pinning energy, the velocity 
tends to increase linearly with the force, V = mF, with 
the mobility m = 1/7, corresponding to the dashed line 
in Fig. [Ha). 

In Fig. Hlb) the disorder intensity effects on the 



velocity-force characteristic can be observed. The crit- 
ical forces for the corresponding disorder realizations for 
each intensity are quoted. Since around the depinning 
transition the velocity strongly depends on the sample 
critical force value, along the present work we will use 
scaled variables for velocity and force. The scaled veloc- 
ity is given by 

which defines a systematics to average over disorder re- 
alizations. Besides, we use as the control parameter the 
scaled force 



which measures the scaled distance to the critical force 
for each disorder realization. These definitions of scaled 
variables are different than the scaled variables used in 
standard critical phenomena. In our case, we are using 
the critical force of each disorder realization in order to 
measure how close the system is to the critical point, in- 
stead of using the disorder averaged value F c . Besides, 
we also incorporate into the definition of the order pa- 
rameter v a non-trivial disorder average when using the 
disorder realization dependent value F c . 

Although a temperature-dependent critical force can 
be considered for studying thermal properties at de- 
pinning (Iol |. we are using here the zero-temperature 
value. From a practical point of view the temperature- 
dependent critical force can be defined as the inflexion 
point of the velocity-force curves. In fact, in the tem- 
perature range we are studying here this temperature- 
dependent critical force does not deviate much from the 
zero-temperature value. Instead of using a temperature- 
dependent critical force, we adhere here to the idea that 
the important quantity given by the disorder potential 
is the zero-temperature critical force and that the small 
temperature data can be interpreted using this quantity. 
Besides, the zero-temperature critical force strictly de- 
pends only on the disorder configuration and therefore 
permits the computation of the average velocity using 
Eq. ([3]) with disorder and temperature averages indepen- 
dently realized. 

The scaled variables, Eqs. (0 and ©, are natural for 
an overdamped particle driven in a periodic potential 
U(x) = R cos(x/\), rydx/dt = —dll(x)/dx + F, where 
one can readily obtain that F c = Ro/X and ^V/Fc k, 
■J (F — F c ) I F c close to F c . They also arise in functional 
renormalization group calculations for center of mass ve- 
locity of an elastic manifold, jV/F c ~ \(F — F C )/F C ]P but 
with 7 an effective friction coefficient [23J . On the other 
hand, from a practical point of view, it was shown that 
the scaled variables defined in Eqs. (|5|) and ©, applied 
to each particular sample, serve to diminish sample-to- 
samplc fluctuations when studying the depinning transi- 
tion of a string [3(1 H3 ■ 
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The inset of Fig. [ljb) shows the same data as in the 
main panel but in scaled form for a single realization. The 
difference between these curves close to the threshold is 
due to the fact that the full function V(F, T) depends 
on the disorder intensity not only through the value of 
the critical force, but one also needs to consider both 
extra disorder-dependent temperature scale and friction 
coefficient. This interesting issue will not be crucial for 
our present study however (see discussion in Sec. IVIII 
below) . 



IV. TEMPERATURE DEPENDENCE OF THE 
VELOCITY AT THE CRITICAL FORCE 

Here we show the finite temperature response of the 
elastic string exactly at the critical force and discuss 
some finite size scaling effects, in particular the crossover 
to single particle dynamics. Figure [2] presents velocity- 
temperature curves for different system sizes. All the 
curves were computed at exactly the sample critical force, 
using the scaled force variable / = 0, and then aver- 
aged over disorder realizations. The disorder intensity is 
i?o = 0.5 and the results are qualitatively similar to those 
reported for R = 1 in Ref. [301. At very high temper- 
atures, T ^> i?o, the system enters the fast flow regime 
and the velocity practically equals the force; therefore the 
reduced velocity (which incorporates the critical force) 
tends to unity. At intermediate temperatures, the ve- 
locity reduces and the curves tend to display the critical 
behavior v ~ T^. This power-law behavior is however 
interrupted by finite-size effects at smaller temperatures, 
when the dynamic characteristic length £ equals the sys- 
tem size L. 

At very small temperatures a crossover to single par- 
ticle dynamics [13, is observed as shown by the 
L = 32 curve. A simple ad-hoc model to rationalize this 
crossover has been given by Duemmer and Krauth [37l ] 
while numerically studying the zero-temperature depin- 
ning transition. Within this model, one can write the 
velocity in the very small temperature regime as 



M 



T +n(T)' 



(7) 



where t±(T) is the temperature-dependent time the in- 
terface spends near the critical configuration and tq is 
the rest of the time needed to cover the transverse spa- 
tial period M of the computational box. In this simple 
model, To is approximated to be temperature indepen- 
dent at very low temperatures. Using the temperature 
dependence of the escape rate for a particle in a ran- 
dom potential f|cj, one can propose that n = aT~ 1//3 . 
In Fig. [2] we show with a dotted line that the very small 
temperature regime for L = 32 is well fitted with Eq. (J7]) . 
For L = 32 and M = 304 as we fitted the T < 10~ 3 

regime using Eq. (JJJ and we found the fitting parameters 
To = 779.5 and a = 3.42. This is a simplified model allow- 
ing to rationalize the crossover to one-particle dynamics 
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FIG. 2. (Color online) Velocity-temperature curves for differ- 
ent systems sizes L, as indicated, while keeping M = 4 L^ dc p , 
with £dc P = 1-25 [13, the depinning roughness exponent. 
The disorder intensity is Ro = 0.5. All the data were com- 
puted at exactly the sample critical force and then averaged 
over disorder realizations. The dashed line corresponds to the 
expected power-law behavior. The dotted line describes the 
crossover to the one-particle regime as discussed in the text. 
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FIG. 3. (Color online) Finite size scaling of the velocity- 
temperature curves for different L values according to Eq. ©. 
The disorder strength is Ro = 0.5 and the transverse size has 
been kept at M — 4L idc p. We show points for T < 0.02 and 
L > 64, since for L — 32 the single particle regime is present 
at small temperatures. The dashed line corresponds to the 
power-law behavior with ip = 0.15. We also used the values 
P = 0.33 [13| and v = 1.33 HJ. 



and should be further tested. 

The finite size effects displayed by the velocity- 
temperature curves in Fig. [2] are not easily accounted for 
by standard finite size scaling arguments. In fact, assum- 
ing finite-size scaling as in standard critical phenomena, 
the velocity should be described by universal functions 
as in 



T 



(8) 



with g(x) ~ 1 for x <C 1 and g(x) ~ for x 3> 1. As 
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mentioned in Ref. [30| strong corrections-to-scaling ef- 
fects are present in these results. In order to show that, 
Fig. [3] presents an attempt to use the standard finite size 
scaling correction scaling, Eq. (JSJ), with the bare data 
in Fig. [2 One can observe strong finite size corrections 
and this can also be observed with other values of Ro . In 
addition, the collapse of the data does not improve signif- 
icantly when using other values of the scaling exponents. 
Despite these strong finite-size effects exhibited by the 
velocity at critical force, the power-law regime charac- 
terized by the thermal rounding exponent ip = 0.15 does 
not suffer from strong finite size effects, as shown in the 
following sections. 



V. STRUCTURE FACTOR ANALYSIS 

In this section we turn to the complementary geometri- 
cal analysis of the structure factor, which contains infor- 
mation on the geometry of the string at different length 
scales. The results presented in this section complements 
those reported in Refs. [3(j,|42[ by including different dis- 
order strengths. 

From the numerical simulations, the steady state struc- 
ture factor is defined as 



L-l 

E u i eiq3 

3=0 



(9) 



where q = 2nn/L, with n = 1, L — 1. One can show 
using dimensional analysis that when the width w of a 
sclf-affinc interface of size L is described through a rough- 
ness exponent £ ; i-e- w ~ L^, then the structure factor 
behaves as S q ~ q~( 1+2 ^) in 1 + 1 dimensions. 

At small length scales, q ^> l/£, the structure factor 
shows the typical roughness regime associated to depin- 
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-(l+2Cde P ) ) wm le at large length scales, 



nmg, 

q <C l/£, fluctuations are dictated by effective thermal 
fluctuations induced by the disorder, i.e. S q ~ q~( 1+2 ^ h \ 
The thermal and depinning roughness ex pon ents are, 



respectively, Cth = 1/2 and £<je 
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In 



the critical region the depinning correlation length is 
given by the velocity as £ ~ v~ v /@. Thus, the depin- 
ning correlation length depends on the temperature only 
through the velocity and in the thermal rounding regime 
£ „ [3(|_ With this information one can write 

for the structure factor that 



(10) 



where the scaling function s(x) ~ x~^ i+2 ^ th ^ for x <C 1 
and s[x) ~ ^-(i+^dep) for x 3> 1. In Ref. [30| we showed 
that the structure factor scales with the previous form 
using L — 1024 and M = L^p for the disorder inten- 
sity Ro — 1. Here, we show in Fig. 0Ja) the tempera- 
ture dependence of the structure factor corresponding to 
Ro = 0.5, L = 1024 and M = 4L ?dt =p. For these parame- 
ters the presented data do not show transverse finite size 
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FIG. 4. (Color online) (Color online) Structure factor and its 
scaling function for Ro = 0.5. The curves in (a) correspond 
to the critical force and different temperatures. The size of 
the sample is given by L = 1024 and M = 4L Cd °p . (b) Scaled 
curves showing the crossover between the depinning regime 
at small length scales (x = qT~^"/P ^> 1) and the thermal 
regime at large length scales (x = qT~^ v ^ <C 1). Together 
with ip = 0.15, the values /3 = 0.33 [13], v = 1.33 and 
Cdcp = 1-25 [33,111] were also used. 



effects [42]. Figure Efb) shows the scaling of the structure 
factor according to Eq. (JTU)) and using tp = 0.15, which 
shows a very satisfactory data collapse. 

In order to reach the steady-state for the same tem- 
peratures of Fig. BJa) but larger disorder intensities it 
is necessary to equilibrate the system for longer times. 
Since this equilibration time scales with the transverse 
system size we can reduce the simulation time by us- 
ing M = L^ dc p for i?o = 5. The resulting data, shown 
in Fig. EJa), presents the small length scale depinning 
regime and the large scale effective thermal regime de- 
scribed above, but also clearly show a larger length scale 
regime where finite transverse size effects are present. In 
this regime the roughness exponent is the one correspond- 
ing to a random-periodic system in the fast-flow regime, 
Cper = 3/2 [42]. Hence, we can detect and discard the 
data corresponding to this random-periodic regime in or- 
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FIG. 5. Structure factor and its scaling function for Ro = 5. 
The curves in (a) correspond to the critical force and different 
temperatures. The size of the sample is given by L = 1024 
and M = Z^ d °p. (b) Scaled curves showing the crossover 
between the depinning regime at small length scales (x = 
qT -^v/f> > ^ and the thermal 

regime at large length scales 
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FIG. 6. (Color online) (a) Short time evolution of the velocity 
at the critical force and for different temperatures, from T = 
0.0001 (lower curve) to T = 0.01 (upper curve), as indicated 
in the key. Data correspond to Ro = 1, L = 1024 and M = 
L^ d °p.(b) Short time scaling of the velocity for t > 20. The 
« 1). The very large length scale random- behavior x~ p/z " for x = tT^ zvlfi < 1 is indicated with a 



periodic fast-flow regime has been discarded (see the text). 
Together with tp = 0.15, the values f3 = 0.33 [13], v = 1.33 gj] 
and £dcp = 1.25 [13, Hit were also used. 



dashed line. Together with tp - 
v = 1.33 EH and z = 3/2 



= 0.15, the values /3 
were also used. 
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der to get a curve that can be scaled again using Eq. (fT0|) 
and the known random-manifold exponents, as shown in 
Fig. [5£b), getting again a very satisfactory collapse. 

Therefore, we have presented here data of the structure 
factor for different disorder intensities which shows that 
the quoted thermal rounding exponent is disorder inde- 
pendent. Furthermore, we have shown how the thermal 
rounding exponent gives the temperature dependence of 
the depinning correlation length, £ ~ t~4>v/0^ from a 
steady-state geometry analysis. 



VI. SHORT-TIME DYNAMICS ANALYSIS 

One possible way to get rid of finite size effects is to 
analyze the short time dynamics. Starting from a given 



non-steady initial condition at fixed force F and temper- 
ature T, the velocity begins to evolve with time until it 
reaches the steady-state value corresponding to the val- 
ues of F and T. This transient dynamics is controlled, 
at short times, by a single growing correlation length, 
£(t), which at longer times saturates to the steady-state 
correlation length above threshold, £ ~ v-v/P. Since 
the transient correlation length grows as £(t) ~ i 1 / 2 , 
with z Ps 3/2 [4lJ the depinning dynamical exponent, 
scaling arguments show that the velocity decreases with 
time as v[t) ~ ^{t)- '" ~ r^ 1 ™ [H[ before saturat- 
ing to the steady-state value above threshold, given by 
v{t ->• oo) ~ f p at T = or by v(t -)■ oo) ~ at / = 0. 

Figure [f^a) shows the time evolution of the velocity 
exactly at the critical force and for different tempera- 
ture values as indicated. The dashed line corresponds to 
the expected short-time critical behavior v[t) ~ t~^l zv . 
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Discarding the very short time regime, t < 20, which con- 
tains information about the microscopic non-universal 
dynamics [4lj . the curves in Fig. EJa) can be recast into 
a universal form using the scaling function 



(a) 



v{t) = T^h [tT^ zv l^ 



(11) 



with h(x) ~ x~^/ zv for i < 1 and h(x) ~ 1 for j; > 1. 
The data collapse shown in Fig. |6jb) uses the previously 
known depinning exponents (3 — 0.33 37], v = 1.33 [4lj . 
z = 3/2 [4l|, together with the thermal rounding ex- 
ponent ip — 0.15. Since the data collapse is good with 
no need of adjustable parameters we can conclude that 
the value of the thermal rounding exponent obtained in 
Ref. is consistent and does not suffer from strong 
finite-size effects. 



VII. VELOCITY SCALING FUNCTION 
AROUND DEPINNING 

In this section we turn to the analysis of the universal 
behavior of the force and temperature dependent veloc- 
ity function, v(f,T). Focusing on testing the robustness 
of the thermal rounding exponent ip, parameter values 
around the critical point given by / = and T — are 
tested. If there were not strong finite size effects, in the 
vicinity of the critical region the velocity should scale as 



HlfT 



-V-//3 



(12) 



with H(x) ~ for i » 1. Figure [TJa) shows velocity- 
force curves for different temperatures and for Rq = 1. 
The numerical data is split into two sets: on the one hand 
data points correspond to given parameters which are 
"inside" the thermal rounding region, and on the other 
hand continuous lines represent data "outside" the ther- 
mal rounding region. The data are outside the critical re- 
gion either because temperature is too high, T > 0.02 in 
the present case, or because the force is far away from the 
critical value, |/| 3> 1 (/ 3> 1 and / -C — 1 corresponding 
to the fast-flow and creep regimes, respectively). In ad- 
dition, to avoid finite-size effects, data points are also 
considered "outside" the critical thermal rounding re- 
gion if they correspond to velocities smaller than the 
crossover at v ~ v m i n to single-particle behavior for each 
size L. Since in the critical region £ ~ v(f,T)~ u /P, as 
shown from the structure factor analysis, we roughly have 
Vmin ~ L~PI V . According to such criteria, the selected 
data is finally presented in the scaled form Eq. (fT2")) in 
Fig. Eth) for / > 0. The dashed line indicates the ex- 
pected asymptotic x^ form, corresponding to the scaling 
function H around the critical region. The collapse into 
a single curve for different T and / confirms numerically 
that the data set used is inside the critical scaling region. 

The scaling function H is not yet universal as it also 
depends on the disorder intensity. Figure [7] shows the 
critical region and the form of H for Rq = 1. In Fig. [5] 



0.5 



1 

G © 


1 1 1 

T= 0.001 


i 




T= 0.002 




□ □ 


T= 0.005 




< < 


r=o.oi 




O O 


r=o.02 






r=o.o5 






^ i , I 






FIG. 7. (Color online) (a) Velocity-force characteristics at 
finite temperatures for Rq = 1. Data shown with points 
(lines) are inside (outside) the thermal rounding region (see 
the text), (b) Scaling of the data for / > and in the ther- 
mal rounding region using the scaling form given in Eq. (|12[) . 
The dashed line indicates the asymptotic expected form 
for x = fT~ i ' /l3 > 1. Here, the value f3 = 0.33 [13 was used 
together with t[> = 0.15. 



we show the velocity scaling function H for different dis- 
order intensities, Rq — 1,2 and 5, for the full force range 
within the thermal rounding region. In Fig. [S] the same 
data is presented in a double logarithmic scale. As can 
be observed all curves display the asymptotic power-law 
form H ~ x@ for / 3> T^l® , but with different prefactors 
for each disorder intensity. 

At this point, in order to properly include the disor- 
der intensity on the scaling velocity function and obtain 
the universal function, a disorder-dependent tempera- 
ture scale T c is needed. Again, for the simple exam- 
ple of the depinning of a particle in a periodic potential 
U(x) = i? cos(x/A), jdx/dt = -dU(x)/dx + F + r)(t), 
with rj{t) a Langevin noise at temperature T, it is easy 
to see, from pure dimensional analysis, that T c ~ Rq, 
and therefore ~/V/F c = h[(F ~ F C )/F C ,T/R ], where 
F c = Rq/X. Such naive approach does not work for the 
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FIG. 8. (Color online) Scaling of velocity-force curves for 
different temperatures and close to the thermal rounding 
regime. Different disorder intensities are shown: Ro = 1 
(blue), Ro = 2 (orange) and -Ro = 5 (black). Symbols rep- 
resent different temperatures as in Fig. EJa). Here, the value 
/3 = 0.33 was used together with if) — 0.15. 



elastic string, as the characteristic energy scale is not 
simply Rq as for the particle, but it arises from the inter- 
play of disorder and elasticity. Although it is not obvious 
that it should work at depinning, one is tempted to use 
the scaled temperature r = T/T c , where T c = U c /k,B 
gives the characteristic energy scale in the creep regime 
at small forces [HI, . This energy scale is however not 
universal and depends on microscopic details of the disor- 
der [l8|,[l{|[43|]. The assessment of the full dependence of 
U c on microscopic parameters is thus not straightforward 
and from a pragmatic point of view one could directly fit 
it from the creep law. As shown with numerical simu- 
lations within the creep regime, at larger temperatures 
than the one used here, U c can be fitted from the creep 
law, but it dependence on R is not trivial [3(|- We do 
not have access to T c with the present numerical results, 
which focused in the force region around the critical de- 
pinning threshold. 

Therefore, we can not incorporate at this stage the 
influence of the disorder intensity in the velocity function. 
In spite of this, the data displayed in Fig. [S] clearly show 
that the velocity can be represented in a scaled form, 
with identical critical exponents, for different disorder 
intensities. More importantly, these data supports the 
disorder independent value ip — 0.15 tested here. 

Finally, it is worth relating our results with the univer- 
sal scaling function proposed by Nattermann, Pokrovsky 
and Vinokur [43| using a phenomenological interpolat- 
ing form for the full force and temperature dependence 
of the velocity of a domain wall in a random medium. 
This form includes the thermal rounding regime around 
F c and the F <C F c creep regime, thus depending also on 
the universal creep exponent /x = 1/4 (for a one dimen- 
sional elastic interface). The proposed functional form 
in Ref. [43| is different below and above the critical force 
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FIG. 9. (Color online) Same data as in Fig. [5] but in double- 
logarithmic representation. The upper and lower curves cor- 
respond to / > and / < 0, respectively. We also show fitting 
curves using a\ +a2X^ and bi exp[— 62 |a?|' s ' , as suggested by 
the scaling functions Eqs. (|15|) and (|16[) . 



and can be written as 



exp 



V = mF- 



-*(!-*)"*<*)' 



1 + 



F ) 



T 



(13) 



for F < F c and 



V = 



mF 



mF 1 - 



1 



21 (£k) 

T \ F I 



Fc 

F 



(14) 



for F > F c . It can be shown that close to the depinning 
region, i.e. above the creep regime, where / = (F — 
F c )/F c <C 1 and t = T/Tj « 1 this phenomenological 
form can be reduced to the scaling form Eq. (fP7|) . with 
H(x) = H~(x) for / < and H(x) = H+{x) for / > 0. 
The corresponding limit functions are 



H- [fr 



-(.-fr 



(15) 
(16) 



Since we do not have the temperature scale T c from 
the creep law, we have directly fitted the data for the 
velocity scaling function to the universal forms sug- 
gested by Eqs. CG2) and (USD- The / < and / > 
ranges have been fitted separately using a± + a^x^ and 
61 exp[— 6 2 |^| /3// '''], respectively, obtaining four fitting pa- 
rameters for each disorder intensity. The results are 
shown in Fig. [5] In all cases the fit is better in the / < 
region. Furthermore, one can observe that the obtained 
curves interpolate badly around / = 0. In fact, enforc- 
ing a\ = b\ makes the fitting considerably worse. We 
therefore conclude that the data can not be satisfacto- 
rily fitted using this phenomenological form, particularly 
above threshold, hence evidently pointing to the need of 
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a more accurate description of the thermal rounding of 
the depinning transition. 

The phenomenological functional forms, Eqs. (fl3l) and 
(H3J), give a potentially important tool which allows to 
directly fit experimental data. This was directly used 
in Ref. [33[ , where the velocity- force characteristic below 
threshold for ultrathin ferromagnetic layers was fitted us- 
ing Eq. (|13|) . By fitting just one experimental curve be- 
low threshold, the value ip = 0.15 ± 0.10 was obtained 
for the thermal rounding exponent. Since several fitting 
parameters were used and due to the large error bar, this 
value can only be compared with our numerical value 
with extreme caution. Anyway, the experimental value 
is consistent with our numerical simulations. 



VIII. SUMMARY 

We have presented extensive numerical simulations to 
test the validity of the thermal rounding exponent of the 
depinning transition. We analyzed the direct scaling of 
the steady-state velocity-force characteristics, the steady- 
state structure factor and the short-time transient dy- 
namics. The existence of a critical (power law) thermal 
rounding of the depinning transition is consistent with all 
our results, together with the existence of a unique di- 
vergent length scale, dependent on temperature and/or 
distance to the critical pinning force, but ultimately con- 
trolled by the velocity as in the zero temperature de- 
pinning transition. The results are all consistent with a 
value of the thermal rounding exponent of ip = 0.15 in 
agreement with our previously reported value [30j |. This 



exponent describes the power-law vanishing of the ve- 
locity with temperature exactly at the critical depinning 
force, V ~ T v , for the universality class of one dimen- 
sional elastic interfaces with short-range elasticity and 
short-range correlations in the disorder. 

Although the value of the thermal rounding exponent 
have been previously obtained with larger system sizes, 
where finite size corrections are still observable, we have 
shown here that this value is also consistent with short- 
time dynamics results which do not suffer from severe fi- 
nite size effects. Besides, ip = 0.15 also describes the scal- 
ing properties of the structure factor for various disorder 
strength values, connecting this value with a geometri- 
cal roughness crossover in the interface. Finally, we have 
shown that it is consistent with a scaling function de- 
scribing the velocity-force characteristics as a function of 
temperature and force. Experimental confirmation of our 
results, directly targeting the thermal rounding regime 
and allowing to test the value of the thermal rounding 
exponent, would be welcome. 
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